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1. Introduction 

Heavy Flavour physics at HERA focuses on as- 
pects related to production dynamics rather than 
weak decays and mixing angles. Heavy quarks 
- like jets - reflect the properties of hard sub- 
processes in ep interactions. Their mass provides 
the natural scale which allows the application of 
perturbative methods in QCD and also testing 
the theory in regions where no other hard scale 
(like high transverse energy) is present. Fur- 
thermore, when probing the structure of matter, 
heavy quarks single out the gluonic content. 

For this talk topics from open charm and 
beauty production have been selected. There was 
no time to cover Onium production (for recent re- 
views, see [0). There is no tau physics included, 
either, although r candidates have been seen at 
HERA, and interesting limits on lepton-flavour 
violating lepto-quarks have been obtained Q. 

At HERA, 820 GeV protons (defining the 
"forward" direction) collide head-on with 27 GeV 
electrons, yielding a centre-of-mass system (CMS) 
energy of ^/s = 300 GeV. The standard textbook 
notation will be used to describe the kinemat- 
ics of deep inelastic scattering (DIS): denotes 
the four-momentum transfer squared, x and y 
are scaling variables related to the momentum 
fraction of the struck parton and to the inelas- 
ticity of the collision, respectively. W denotes 
the photon-proton CMS energy. The kinematic 
regime where the exchanged photon becomes qua- 
si-real (Q^ 0) is called photo-production. Also 
in this region W can be large and attains values 
an order of magnitude higher than in fixed-target 
experiments. 



2. Charm Production 

In most of the studies presented here, open charm 
is detected in the "golden" decay channel 



D°TT+ followed by D° K-tt+} ZEUS also 
uses the channel £)*+ (D" i^^vr^Tr+Tr")??", 
and this year for the first time they have shown 
resuhs using ^ (0 ^ K+K-)tt+ ||. 

2.1 Deep inelastic scattering 

In QCD, c and b production in ep collisions pro- 
ceed mainly via the boson gluon fusion diagram 
shown in figure |l|. 
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Figure 1: Boson gluon fusion. 

Therefore, charm production has already in 
the past offered a way to determine the gluon 
density in the proton The process has been 
calculated in Next-to-Leading order (NLO) QCD 
in the so-called Three Flavour AIS scheme for 
photo-production § and DIS §. In the HVQDIS 
program, the charm hadronization into D* me- 
sons is modeled using a Peterson fragmentation 



^The charge conjugate is always implicitly included. 
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Figure 2: D* cross sections in DIS, compared with NLO QCD, using for fragmentation a Peterson type 
(open band) or a LO MC model (shaded band) . The widths of the bands correspond to a charm mass variation 
between 1.3 and 1.5 GeV. 



function [Q with the parameter Cc = 0.035 as 
determined in e"'"e~ annihilation 

Differential D* cross sections in the exper- 
imentally accessible DIS range have been mea- 
sured by HI § and ZEUS The ZEUS data 
are shown in figure^ and compared to the HVQDIS 
calculation. (The open points in (d) and (e) are 
results from the KAtt channel.) A Leading-Order 
Monte Carlo fragmentation approach based on 
the JETSET program also has been used by ZEUS. 
As can be seen from the figure, the pseudo-ra- 
pidity and the xd distribution^ are quite sen- 
sitive to details of the fragmentation modeling, 
whereas the distributions of other kinematic vari- 
ables, like X and are barely affected. The mi- 
gration of the mesons towards positive pseudo- 
rapidities in the Monte Carlo, due to color inter- 
actions between the charm quark and the proton 
remnant, has been called the "beam drag" ef- 
fect 0. 

In general, there is good agreement with the 

^xo = 2p*{D*)/W is the fractional momentum of 
the meson in the ■y'p rest frame and strongly correlated 
with r]. 



NLO calculation. This gives the justification to 
the extrapolation to the full phase space, which 
is needed to extract , the charm contribution 
to the proton structure function. 

In the Quark-Parton Model, the structure 
function i*2 is given by the charge-weighted sum 
of the (anti-)quark densities, F2 = J2i ^liQii^) + 
qi{x)) and depends only on x. The presence of 
gluons introduces "scaling violations" , i.e. a de- 
pendence of the structure function on . The 
gluons are also responsible for the observed steep 
rise of F2 towards low values of x, where the 
quark sea is being probed. 

The relative charm contribution to inclusive 
DIS, quoted as F^jFi, is not constant, but rises 
with and towards low x, as can be seen in 
figure 1^. Since the process under study is gluon- 
induced, F^ reflects the gluon content of the pro- 
ton even more pronouncedly than F2. The most 
salient feature is that in the HERA regime the 
charm contribution is very large, 20 - 30%, mak- 
ing the theoretical description of charm produc- 
tion an essential ingredient to the understanding 
of proton structure. 
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Figure 3: Relative charm contribution to inclusive 
DIS. (Extrapolation uncertainties of ^ 20% are not 
shown.) 

HI has performed a direct measurement of 
the gluon density in the proton from charm pro- 
duction in DIS and in photo-production. The ki- 
nematics of the reconstructed D* meson and the 
scattered electron have been used to infer the in- 
coming gluon momentum. The NLO programs 
have been used in order to correct for higher or- 
der processes and fragmentation by means of an 
unfolding procedure. The results, shown in fig- 
ure |[ agree well with each other and with the 
gluon distribution determined indirectly through 
a QCD analysis of the scaling violations of inclu- 
sive structure function data. This demonstrates 
that the gluon density is universal, and indicates 
the success of the NLO QCD description. 

2.2 Photo-production 

In the case of photo-production of open charm, 
additional, so-called "resolved" photon contribu- 
tions have to be taken into account. The photon 
may fluctuate into a hadronic state, and a par- 
ton from this state interacts with a parton from 
the proton, e.g. charm may be formed via gluon- 
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Figure 4: Gluon distribution. 

gluon fusion. Since only a fraction of the pho- 
ton's energy is available in the hard sub-process, 
in comparison with "direct" photon interactions, 
the outgoing charm is found at lower transverse 
momenta and at more forward rapidities. Re- 
solved photo-production is therefore suppressed 
by experimental cuts, and in fact extra care was 
taken in the above-mentioned gluon analysis to 
minimize uncertainties related to such a priori 
unknown contributions. 

There are two approaches in NLO QCD to 
describe charm photo-production. In the "mas- 
sive" scheme Q , charm is produced only dynami- 
cally in the final state (like in the DIS case above). 
In the "massless" scheme charm also plays 
an active role in the initial state as a parton of 
the proton and of the photon. While the first ap- 
proach should be adequate near threshold where 
the transverse momenta are comparable to the 
charm quark mass, the second has been devel- 
oped for the region of higher p± . 

D* photoproduction cross sections as recently 
measured by ZEUS § in the region 80 < < 
120 GeV are compared to NLO QCD calculations 
in the massive scheme in figure ^. Fair agree- 
ment can be seen in the low p±^ range, with how- 
ever somewhat "stretched" parameter values. At 
higher p± the data tend to be above the predic- 
tion, in particular in the forward region. As in 
the DIS case, predictions in this region are seen 
to be sensitive to details of how the fragmenta- 
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Figure 5: D* photoproduction cross section, com- 
pared with NLO QCD in the "massive" scheme. 
(Dashed: Peterson, dotted: PYTHIA fragmenta- 
tion). 



tion process is being modeled. The sanie conclu- 
sions can be drawn from the ZEUS data obtained 
with D* and Dg mesons at higher W HI, 
albeit with somewhat larger measurement errors, 
do not observe a discrepancy with respect to the 
"massive" calculations Q. 

In the "massless" calculation, there is a large 
"resolved" contribution which is mostly due to 
diagrams where a charm quark from the photon 
interacts with a parton from the proton; an ex- 



ample is shown in figure |g. (Note that the assign- 
ment of "direct" and "resolved" is not unambigu- 
ous in NLO.) In this picture the cross sections are 
thus sensitive to the charm density in the photon, 
as indicated in figure ^ where different photon 
PDFs have been used in the calculation. HI |^ 
and ZEUS data (shown here) have been com- 
pared with NLO QCD in the "massless" scheme. 
As in the "massive" case, the agreement is not 
good everywhere. 

2.3 DifFractive production 

Pursuing these concepts further, charm produc- 
tion may also shed light on the dynamics of diffrac- 
tive scattering. A fraction of about 10 % of the 
DIS events have a "rapidity gap": in contrast 
to the general case the region surrounding the 
outgoing proton beam is void of any particles. 
This distinctively "diffractive" topology leads to 
an interpretation of the events in terms of the 
exchange of a color-less object, which may be 
identified with the Pomcron (IP) in the frame- 
work of hadron-hadron interaction phenomenol- 
ogy. HERA offers the possibility to study these 
phenomena in a perturbative regime and to in- 
vestigate the partonic structure of the exchange, 
where charm production again singles out the 
gluonic component. 

Two different model predictions are consid- 
ered here (figure]^). In the "resolved Pomeron" 
model Q, the structure of the exchanged object 
is described in terms of quark and gluon densi- 
ties which have been obtained from a fit to inclu- 
sive diffractive DIS data . In the "two gluon" 
model |l^, the color-less exchange is realized in 
terms of two hard gluons, and cross sections are 
sensitive to the gluon density in the proton. The 
two approaches lead to remarkably different kine- 
matic distributions. Both HI and ZEUS |l|] 
have presented first measurements of differential 
cross sections this year. 

zp is an observable correlated with the frac- 
tion of momentum of the exchanged object which 
is carried by the parton interacting with the cc 
pair. In the "2g" model, zp = 1 would hold if the 
partons could be directly observed. The data Q 
in figure ^ reveal that a sizeable fraction of charm 
is produced at low z^*"* ( "resolved Pomeron inter- 
action"). The "resolved IP" approach however 
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Figure 8: D* cross section in diffractive DIS. 



P P' 

Figure 7: Diagrams of diffractive D* production: 
"resolved Pomeron" (upper) and "2 gluon" model. 

fails in reproducing the normalization of the HI 
data. In the "two gluon" model, on the other 
hand, a lack of higher order, low zp contribu- 
tions is apparent. More general, models, where 
the hadronic system X predominantly consists 
of the cc system alone, are disfavored. The data 
still have large, mostly statistical, uncertainties, 
but they already provide valuable and rather dis- 
tinctive indications for further refinement of the 
theoretical description. 



3. Beauty Production 

Finding beauty at HERA is not straightforward. 
Theory predicts total cross section ratios of about 

^uds ■ f-^charm ■ f^bcauty ^ 

2000 : 200 : 1 ||. The 

measurements to date rely on the well established 
signature of semileptonic decays of b hadrons in 
jets. Due to the higher b mass, the leptons in b 
events have a higher transverse momentum p!^'' 
with respect to the jet direction which approxi- 
mates the flight direction of the decaying hadron 
(figure The b cross section can thus be ex- 
tracted by means of a fit to the pJ^^ distribu- 
tion. 

The preliminary HI result released in 1998 
has meanwhile been published The visible 

cross section, quoted for the range < 1 GeV^, 
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Figure 9: b production signature and spectrum. 



0.1 < y < 0.8, > 2GeV, 35° < 6"^ < 130° of 



all' = 0.93 ± 0.081^- nb is found to be surpris- 
ingly high: the Monte Carlo prediction based on 
Leading Order QCD calculation is 4 to 5 times 
lower (c^^'* = 0.191 nb). NLO corrections are 
however non-negligible, as in the charm photo- 
production case. 

In figure |l^ the result, extrapolated to the 
full phase space by means of a NLO program 
incorporating a Peterson type model of fragmen- 
tation, is confronted with the NLO prediction. 
Theory at NLO still undershoots the experimen- 
tal result by about a factor of 2. This is in 
marked contrast to the expectation that predic- 
tions for beauty should be more reliable than for 
charm, due to the higher scale set by the b mass. 
On the other hand, the NLO approximation also 
has difficulties reproducing the normalization of 
open b cross section measurements at the Teva- 
tron [|0). 

ZEUS has determined the beauty cross sec- 
tion in three independent analyses, using elec- 
trons, muons in the central region and muons 
in the forward detector region The p^'' 

distribution, obtained in the e channel, is dis- 
played in figure ||. Overlayed is a Monte Carlo 
prediction obtained with the HERWIG genera- 
tor, which already describes the shape of the 



data reasonably well. It should be noted that 
in HERWIG about half of the c and b photo- 
production is due to resolved processes with a 
heavy quark in the initial state, mostly as an 
active parton in the photon. A fit of the dif- 
ferent components yields a b fraction of 20%, 
which translates into a visible cross section of 
(Tvzs (e+ 2 jets + e- +X) ^39± UtH pb for 
events with 2 jets {E^^ > 7, 6 GeV, \r]\ < 2.4) in 
the kinematic range < IGeV^, 0.2< y <0.8, 
e: Pt > l.eCeV, \ri\ < 1.6. 

Consistent results are obtained in the muon 
channels, with the same jet requirements and for 
the same kinematic range. For central muons 
(-1.75 < ?7 <1.3, p >3 GeV) ay^s{e+ 2 jets + 
fi + X) = 36.4 ± 5.2tl°i^ pb is quoted, and for 
forward muons (1.4< rj <2.4, p± >1.5 GeV) 
<Jvts{e+ 2 jets + i^i + X)^ 20.5 ± G-btfj pb. 

Using the Monte Carlo to convert the mea- 
sured hadronic cross section to a partonic cross 
section allows a comparison with NLO QCD cal- 
culations at parton level. This is done for the 



ZEUS result in the e channel in figure |10[ Here, 
the result is found to be even a factor of 4 above 
theory, with somewhat larger errors. 

The long lifetimes of c and b hadrons which 
can be measured with microvertex detectors pro- 
vide an independent signature. With its Central 
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Figure 10: b production cross section, compared with NLO QCD. 



Silicon Tracker (CST), HI has for the first time 
at HERA attempted to exploit these features. 
An event with two jets and a penetrating track 
identified as a muon is shown on the left-hand 
side of figure H]. The muon has a relatively high 
transverse momentum of 3 GeV/ c with respect to 
the nearest jet, which is unlikely for muons from 
charm decays or mis-identified hadrons but is ex- 
pected for decays of b flavoured hadrons. In the 
magnified view perpendicular to the beam (on 
the right-hand side) the tracks measured in the 
CST are represented as bands with widths cor- 
responding to their ztlcr measurement precision. 
The resolution provided by the CST reveals that 
the muon track originates from a well separated 
secondary vertex. This is now being integrated 
into ongoing physics analyses. 

4. Summary and Outlook 

Heavy quarks have become an increasingly in- 
teresting part of the physics field opened up by 
HERA. With the statistical power now reached in 
the data, charm quarks are a direct probe of the 
gluonic component of hadronic structure. In the 
case of the proton, the fact that about one quar- 
ter of the ep interactions in the HERA regime 
result in final states with charm makes it evi- 
dent that understanding the charm contribution 
to the structure function F2 is a conditio sine qua 
non. The agreement of the gluon density directly 



determined from charm with the result from the 
scaling violation analysis lends strong support to 
the underlying QCD picture. 

New light has also been shed on the struc- 
ture of the photon, where a charm content in the 
"resolved" photon may provide a viable route to 
understanding the HERA data. Finally, probing 
the structure of diffractive exchange with charm, 
differential cross section data have started to dis- 
criminate between different theoretical concepts. 

Beauty production at HERA has now been 
measured by both HI and ZEUS, and found to 
occur at a rate that is presenting NLO QCD with 
a challenge. Although experimental errors are 
still considerable, a re-evaluation of theoretical 
uncertainties might be appropriate. 

Top pair production is kinematically exclu- 
ded at HERA energies. It has however been spec- 
ulated recently that the events with isolated 
leptons and large missing p±, of which the HI 
collaboration observes more than expected 
(and some with atypical kinematics), might be 
due to single top production mediated by a new, 
flavour-changing effective interaction. The top 
quark would decay via t — > Wb, W fiiy, giv- 
ing rise to the observed signature. In fact, some 
of the observed outstanding events are kinemat- 
ically not inconsistent with such a hypothesis. 
Lifetime-based tagging techniques will in the fu- 
ture be used to investigate such exciting possi- 
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Figure 11: A 6 candidate in tlie HI detector; right: vertex region. 



bilities further. One more reason why HI and 
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